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Abstract Near a Feshbach resonance, the two-body scattering length can assume any value. When 
it approaches zero, the next-order term given by the effective range is known to diverge. We consider 
the question of whether this divergence (and the vanishing of the scattering length) is accompanied 
by an anomalous solution of the three-boson Schrodingcr equation similar to the one found at infinite 
scattering length by Efimov. Within a simple zero-range model, we find no such solutions, and conclude 
that higher-order terms do not support Efimov physics. 



1 Introduction 

In a three-body system of identical bosons, an infinite sequence of bound states appears at the threshold 
for binding of the two-body subsystems as shown originally by Efimov in the context of nuclear physics 
QUE]. However, such states unfortunately do not exist in nuclear system [5J. Efimov states have been a 
focus of the ultracold gas community since the experimental observations of signatures of such states in 
a dilute gas of Cesium atoms [HG3IB] ■ This was followed by studies in three-component 6 Li [T1151I WT011T1"] , 
bosonic 7 Li [T211T3"] and 41 K [14], and heteronuclear mixtures of 41 K and 87 Rb [15]. By sweep a Feshbach 
resonance in 7 Li over several orders of magnitude [TO] a recent study [T7] was able to identify eleven 
distinct features in the recombination rates. Some of these are interpreted as four-body resonances 
that have been theoretically predicted as associated to the Efimov trimers [TF1I19U20"] . Evidence for 
potential four-body resonances was also found in [T4"l[2l"j . The universal theory of four-body states has, 
however, been questioned [22] . and it is unlikely that universality persists to higher particle numbers 

The characteristic geometric energy and size scaling of Efimov states near a resonance has been 
discussed within various models by many people (see |24] and [25] for comprehensive reviews) . In cold 
gases one can use Feshbach resonances to tunc interactions 26,27,28,29 and thus reach resonance 
conditions in a controlled way. This gives the celebrated scaling factor e _7r ™/ s ° « 22.7 with sq = 
1.00624. The energy of the nth state is then simply related to the first one by E(n)/ E(0) = 22.7~ 2n 
and the root-mean-square size by r rms (n) /r rms (0) = 22. 7~". This result is found within the universal 
theory where the scattering length is much larger than the range of the interatomic potential. However, 
evidence for systematic deviations from the universal scaling has been found [14] . and therefore it is 
necessary to consider correction to the universal picture in detail. The experiment reported in Ref. [14] 
was found to agree qualitatively with the corrections found in [501131 j . 

In this paper we study the corrections to Efimov physics as the width of the Feshbach resonance 
becomes smaller. In this case the effective-range is sizable and must be taken into account. The zero- 
range model is attractive due to its simplicity and range corrections have been discussed in Rcfs. [311 
32,33,34,35,36 . Here wc use a two-channel model to describe the Feshbach resonance [57] . This kind 
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of model has been considered by several authors for both broad and narrow resonances in few- and 
many-body contexts [MHMllira^ A two-channel model gives 

the energy dependence of the interaction in terms of the background scattering length and the width of 
the resonance. This is essential since we will be interested in the behavior of the corrections both around 
the resonance and around the point where the scattering length goes to zero. We recently considered 
the effects of Feshbach resonances around zero-crossing for a macroscopic Bose-Einstein condensate 
[541155) and found considerable changes in critical particle numbers for very narrow resonances. 

Here we consider the interesting question of whether higher order corrections tied to the effective- 
range term can produce anomalous solutions to the three-boson problem as seen around the resonance. 
While previous studies using finite-range models [35] and full numerical three-body calculations [35] 
have considered similar situations, we simplify the discussion considerably by using zero-range model 
with a two-channel description of the scattering process. This recovers known results about both broad 
and narrow Feshbach resonances. In addition, we investigate the solutions around the point where the 
scattering length goes to zero. Here we find no evidence for anomalous solutions that could produce a 
spectrum similar to the one found by Efimov, regardless of whether the resonance is broad or narrow. 

2 Feshbach Model 

Since we are interested in finite-range corrections within the zero-range approximation, we need to 
consider multi-channel models. In the case of s-waves, two-channel models have been derived using 
various methods, including low-energy effective theory [37], multi-channel quantum-defect theory [5j)l 
[57], and resonance models [3j5]. The on-shell open-open channel s-wave T- matrix as a function of 
magnetic field strength, B, can be written 



where m is the atomic mass, k is the center-of-mass momentum, at> g is the background s-wave scattering 
length, Bo is the position of the resonance, AB is the resonance width, and Afi is the difference in 
magnetic moment of the open and closed channels. We can compare this to the usual expression for 
the T-matrix in vacuum 




(1) 



k cot 5(k) — ik 

to get an expression for kcotS(k) as a function of B. This yields 



T(k) 



m 



(2) 





If we neglect the /c-dependent term we recover the usual formula 



a(B) 



(4) 




Keeping the k term and going to the resonance where B = Bq wc find 



k cot S(k) 



1 2ft 2 

2 mA^ABa^g 



(6) 



1 + 2 



2rnA[iABabg 



a bg k 2 



If we define r e o = —2h 2 /mab g A^,AB and make a low momentum expansion we see that this is consistent 
with an effective range expansion where the constant term is absent as the a{B) diverges at B = Bq. 
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Notice that a bg ApAB > for all Feshbach resonances [29] and therefore r e o < always. Similarly, at 
zero-crossing where B = Bq + AB, we find 



kcot5(k) 



1 

a bg 



a lg r eO 



(7) 



We thus see that no meaningful effective range expansion is possible. However, the T-matrix stays 
finite and one can define an effective potential at zero-crossing that involves a bg and r e o [58] . If we 
introduce the function f(B) = (B — Bq)/AB we can write the general form 



k cot 6(k) = 

a bg 



f{B)-\a bg \r e0 \k 2 



I- f{B) + \a bg \r e0 \k 2 



OS) 



where r e o < was used. We recover the expression above at resonance where /(Bq) = and zero- 
crossing where /(-Bo + AB) = 1. An expansion of the phase-shift to second order gives an expression 
for the magnetic field dependent effective range of the form [SB] 



1 - 



a bg 



-i 2 



a(B) 



(9) 



which clearly shows the divergent nature of r e (B) when a(B) — > 0. Higher-order corrections do, how- 
ever, depend on the combination a 2 r e which has a finite limit |541l55j . 



3 Effective Hyperspherical Potential 

We will employ the hyperspherical formalism |24tl25[[59] , where Efimov states are found from an eigen- 
value equation. In the universal limit where only the scattering length is important we have 



-i/cos(i/f ) + ^= sin(i/f ) ^2 p 
sin(^^) a 



(10) 



Here p is the hyperspherical radius defined for three identical particles by p 2 = \{r 2 2 + r f 3 -fr^) where 
Tij is the distance between particles i and j |24j . For later convenience, we define the function g(v) to 
be the left-hand side of Eq. (flU)) , and also 

-scosh(s^) + -§=sinh(s|) 

g(s) i Tr ^ (ii) 

smh(s-2) 

which is obtained from the left-hand side of Eq. (flU)) by the substitution v = is. The solution of this 
equation for each value of p gives the effective hyperradial potential 

where Q{p) is a diagonal coupling term in the hyperspherical formalism [24]. At the resonance where 
\a\ = oo, we find a constant solution v = iso where sq = 1.00624 (which gives the scaling parameter 
discussed in the introduction). In this case the hyperradial equation for large distances reduces to 

92 4 + 1 1 k 2 ) Hp) = o, (13) 



2 



dp 2 p 

where —h 2 n/2m = E is the energy of the three-body state. Here we have neglected Q{p) at large 
distance; i.e. ro«/)« \a\ [24], where r is a small distance cut-off (sometimes refered to as a three- 
body parameter). In the rest of this paper we will be interested in the asymptotic region ro p. When 
there is a two-body bound in the system (as we will assume below), the leading order in Q(p) behaves 
as —1/Ap 2 at short distance, while in the absence of such two-body bound state, Q(p) falls off faster 
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than 1/p 2 [321160] . In any case, we are interested in the occurence of the Efimov effect which hinges on 
the coefficient of the 1/p 2 in the region ro <C p, so we neglect the contribution from Q(p) from now on. 

Introducing scaled variables fo = /o/ ^fp and p = up, Eq. ()13[) becomes a Bessel equation. The 
solution is ^fpKi SQ {np) where K is the modified Bessel function of the second kind. As the energy goes 
to zero we recover the Efimov spectrum which arises because Ki SQ {p) is log-periodic at low energies 
and thus allows for infinitely many nodes |24] . Below we will be concerned with changes in the Efimov 
scaling as we vary the parameters of the two-channel Fcshbach model. 

In the above we have neglected couplings to higher partial waves which are generally small, and we 
assume that the well-known Thomas collapse is remedied by a boundary condition at small hyperradius; 
/o( r o) = 0. The latter condition in fact determines the position of the first Efimov state and calibrates 
the spectrum. It is not clear whether the actual value of ro can be determined in the universal theory or 
not. Recent measurements seem to suggest that ro is simply proportional to ?* v dw [61] but the reasons 
are not entirely understood (see Ref. |62| for a simple argument as to why this should be true). 

The next step is to replace the scattering length in Eq. (fTU]) by the full phase shift. In order to do this 
we have to use the correspondence between k and the hyperradius p which is k — > k p = v{p) /(\^2p) [321 
[33] . We can then put k p cot 6(k p ) on the right-hand side of Eq. (flTJj) . where S(k p ) is the phase-shift of the 
scattering of a pure two-particle system evaluated at momentum k p . This approach was already used 
in |32U34| . However, as discussed in Ref. [35] . it is in fact an approximation. The potential one should 
use instead is the modified two-body term in the hyperspherical formalism V p (r) = V (\/2/9 sin(-^^)) . 

As was shown in |35j . for an effective range expansion this gives corrections to the effective potential 
as compared to using the standard two-body phase shift and the result in fact deviates from the 
predictions of effective field theory [3T|. The corrections to the effective potential are of order 1/p 3 
and in the large distance regime that we are currently interested in we will neglect this and use the 
simple approximation with S(k p ). We expect that an analysis along the lines of [35] can be made for 
the two-channel model as well and it would be interesting to see the quantititative nature of this in 
the future. 

Replacing — 1 /a in Eq.E0]by Eq. © with k p = v/{V2p), we find 



sO) = — 

a bg 



f(B)-\a bg \r e0 \^ 
l-f(B) + \a bg \r e0 \ 



(14) 



We will now proceed to solve this equation for v(p) and compare to the universal case for different 
values of B. The universal Efimov scaling of energies and radius emerges in the limit of constant 
v = isQ. When the potential deviates from the constant value it is relevant to ask over what region (if 
any) v remains essentially flat and what the numerical value is that gives the scaling. The number of 
bound states is then roughly proportional to ln(pi/po) where po and p\ are then boundaries of the flat 
region [TUS]. 



4 Analysis 

We will consider the case of a bg > only, the opposite case is similar. This implies that there is a two- 
body bound state in the scattering channel. This is reflected in the hyperradial potential for p < a bg 
as discussed in Ref. (32] ■ Since we are interested in three-body physics only, we will consider the region 
p > a bg only. Having introduced the function f(B) above, we can split the analysis of Eq. (|14[) into 
two cases, i) near resonance, f(B) = 0, and ii) near zero-crossing, f(B) = 1. We proceed to analyse 
these in turn. 

4.1 Near resonance 

Here we have f(B) — > 0. The right-hand side of Eq. (|T4"| becomes 



V2p 


4 abg \ 


7) 


2 " 


a bg 


4 a bg 


(* 


rJ 



(15) 
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First consider the case at g 3> |r e o|- Then for any hnitc p, Eq. (|15[) will go to zero. This will also be true 
for the case where p 3> a,b g . Either way we have g(v) = which is known to have the anomalous solution 
v = isg as discussed above. We thus have Efimov trimers in the region bounded by ab g <C p <C \a\ as 
usual. For shorter distances, we have two-body physics at play. In the case a^g ~ |T e o|, the conclusion 
is the same. 

Next consider the case where |r e o| 3> a& 9 which is the case for narrow Feshbach resonances [36,29 . 
For any finite p < \r e o\, we find instead of Eq. (|T5|) . that the right-hand side of Eq. (|T4"|) will behave 
like —\plpja-bg. Looking at Eq. (fTTj) . we see that g(s) becomes a linear function of s for s > 1. A 
solution of the form v{p) = aip/atg can therefore be found in the region db g <C p < |r e o|. However, the 
effective potential in Eq. (|12p then gets a constant contribution and this will not produce an Efimov-like 
spectrum since the requirement is that v 2 be constant and negative over an extended region. 

The final region to consider is p 3> | r* e o | - Here we see from Eq. (|15[) that we recover g(v) = and 
the condition for the Efimov effect is fulfilled. The region is [r e o | "C p < \a\ as previously found [32,35, 
136] . The zero-range model is thus seen to reproduce this intuitively clear result in a simple and elegant 
manner. 



4.2 Near Zero-crossing 



Next we consider the case of zero-crossing where a — > 0, f(B) — > 1, and the right-hand side of Eq. (fl"4"| 
becomes 



\a b gj 



1 



1 keol t ,2 I Obg 

4 «i 3 \ P 

keol t .-2 
a-bg 



2 1 



(16) 



First consider the case where a^g 3> |r e o|. In this case, Eq. (fTB| blows up because of the \r e0 \/a,b g 
term in the denominator. Neglecting the second term in the numerator (which is small since we are 
still assuming that p 3> db g ), we can multiply Eq. (|14[) by v 2 and in turn repeat the argument above 
about the behavior of g(s) for large values of s to obtain a solution that goes like v oc ifip/abg where /3 
contains a factor (abg/jreol) 1 / 3 which we assumed to be large. Again we are lead to the conclusion that 
the effective hyperradial potential in Eq. (fl"2"j) will only receive a constant contribution in the region 
|r e o| dbg </)< \f{B) — 1| • The upper bound is given by how close we are to the point where a 
vanishes, similar to the case near resonance where a is the upper bound. At zero-crossing, the interval 
becomes unbounded from above, and no Efimov-like potential will emerge anywhere for resonances 
where db g \r e o\. The same is true for ab g ~ |r e o| by similar analysis. 

Now consider instead the opposite limit |r e o| 3> a& g . For p <C l^eo), Eq. flU]) becomes simply 
— splpjdbg. The solution in this case is again v oc ijp/cibg and nothing interesting happens. The last 
regime where one can still hope for a non-trivial solution is when p 3> \r e o\, which is where Efimov 
trimers appear near resonance. However, in similar fashion to the case where a,b g 3> |^ e o|j we find that 
v cx iSp/abg with S containing a factor (abg/lreol) 1 / 3 which is now a small number. This does not, 
however, change the outcome and no Efimov-like solution is found. 

We have done numerical root finding on Eq. (fT4")l to confirm the analysis above. We have tested 
the numerics and made sure that it reproduces the usual Efimov solution, sq, and also that it finds the 
free solutions |24| . In the case of zero-crossing, we have not found any constant and purely imaginary 
solution for v. Thus, numerically we also do not see any signs that Efimov-like features are present 
around zero scattering length. 



5 Summary 

We have considered the occurence of a geometric spectrum in the three-boson system around Feshbach 
resonances in ultracold gases when effective-range corrections are included. In particular, we have 
explored the possibility of having an Efimov-like effect around the zero-crossing where the scattering 
length goes to zero and higher order effective-range terms must be considered. This is done within a 
simple zero-range model that allows one to obtain a closed equation for the coefficient of the effective 
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potential in the hyperspherical approach. The effective potential in turn determines whether there can 
be universal three-body states with geometric scaling properties. 

In the analysis of the closed equation, we recover the conclusions of previous studies, including 
those obtained through numerical approaches. Close to the resonance, the region of space where uni- 
versal trimers are supported is bounded from below by either the background scattering length or the 
background effective range, whichever is larger, while it is bounded from above by the value of the 
scattering length (which goes to infinite on resonance). 

Around the point where the scattering length crosses zero, the zero-range model predicts that no 
region can be found where the effective potential can support Efimov trimer states. This is true for 
both broad and narrow Fcshbach resonances. Thus even though the magnetic field dependent effective 
range of a two-channel model diverges around zero-crossing, this does not seem to support a geometric 
spectrum of Efimov trimers. More generally, this implies that higher-order terms do not add attraction 
in the three-body channel of a many-body system. This is consistent with the fact that the instabilities 
observed in condensates with higher-order interaction terms in Rcfs. [54 ( 55 originate from two-body 
effects. 
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